First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.     

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.     

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost , a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.     

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.     

辽东丁香完整叶绿体基因组的结构与特征

为获得辽东丁香（Syringa villosa subsp. wolfii）叶绿体全基因组的基本特征,采用高通量测序技术分析了其叶绿体基因组序列信息,并讨论其系统演化位置。结果表明：（1）辽东丁香叶绿体基因组全长156 517 bp,具有典型的四分体结构;具有131个功能基因,包括36个tRNA基因、8个rRNA基因和87个蛋白质编码基因。（2）该叶绿体基因组蛋白编码区的总密码子偏好性（RSCU）分析显示,RSCU值>1的密码子有31个,其中以A/U碱基结尾的有21个;RSCU值<1的密码子有34个,其中以G/C碱基结尾的密码子有22个。（3）在辽东丁香的叶绿体基因组中,检测出334个散在重复序列,包括170个正向重复序列和164个回文重复序列;检测到227个SSR位点,其中226个位点成功设计出PCR引物。（4）最大似然法构建系统进化树分析显示,辽东丁香与云南丁香 （S. yunnanensis）亲缘关系最近。本研究通过对辽东丁香叶绿体基因组重复序列、IR边界、系统发育等进行分析,为辽东丁香后续的分子标记开发、系统发育分析、物种资源鉴定评价、DNA条形码开发等提供参考。     

Mitochondrial genome phylogeny reveals the deep-time origin of Gomphomastacinae (Orthoptera: Eumastacidae) and its alpine genera in China

Gomphomastacinae is a grasshopper subfamily in Eumastacidae, with a morphology and distribution distinct from other subfamilies. The alpine genera of Gomphomastacinae that inhabit the Qinghai–Tibet Plateau in China show unique characteristics adapted to high-altitude life. However, their phylogenetic position and biogeographic history remain controversial. Thus, to determine the diversification history of these alpine genera and the origin of the subfamily, we obtained mitochondrial genome sequences from all seven Gomphomastacinae genera distributed in China. The reconstructed phylogeny was well supported and confirmed the phylogenetic position of Gomphomastacinae within Eumastacidae. Time calibration revealed a deep-time origin of the subfamily dating back to the Cretaceous period, and the diversification among alpine genera was also an ancient pre-Miocene event (30–50 Ma). Based on phylogeny and time estimates, the most likely biogeographic scenario is that Gomphomastacinae originated from an ancestral lineage that lived in East Gondwana and dispersed to Central and Western Asia through India. Subsequently, the alpine genera likely diverged along with the uplift of the Qinghai–Tibet Plateau and survived drastic climate change by in situ adaptation to high-altitude dwellings.     

Comparative analysis of 343 plastid genomes of Solanum section Petota: Insights into potato diversity,phylogeny, and species discrimination

Common potato (Solanum tuberosum L.) and its wild relatives belong to Solanum section Petota. This section's phylogeny and species delimitation are complicated due to various ploidy levels, high heterozygosity, and frequent interspecific hybridization. Compared to the nuclear genome, the plastid genome is more conserved, has a haploid nature, and has a lower nucleotide substitution rate, providing informative alternative insights into the phylogenetic study of section Petota. Here, we analyzed 343 potato plastid genomes from 53 wild and four cultivated species. The diversity of sequences and genomes was comprehensively analyzed. A total of 24 species were placed in a phylogenetic tree based on genomic data for the first time. Overall, our results not only confirmed most existing clades and species boundaries inferred by nuclear evidence but also provided some distinctive species clade belonging and the maternally inherited evidence supporting the hybrid origin of some species. Furthermore, the divergence times between the major potato clades were estimated. In addition, the species discriminatory power of universal barcodes, nuclear ribosomal DNA, and whole and partial plastid genomes and their combinations were thoroughly evaluated; the plastid genome performed best but had limited discriminatory power for all survey species (40%). Overall, our study provided not only new insights into phylogeny and DNA barcoding of potato but also provided valuable genetic data resources for further systematical research of Petota.     

Genomic analyses reveal natural selection on reproduction related genes between two closely related Populus (Salicaceae) species

Identifying the factors that cause reproductive isolation and their relative importance in species divergence is crucial to our understanding of speciation processes. In most species, natural selection is commonly considered to play a large role in driving speciation. Based on whole genome re-sequencing data from 27 Populus alba and 28 Populus adenopoda individuals, we explored the factors related to reproductive isolation of these two closely related species. The results showed that the two species diverged ~5–10 million years ago (Ma), when the Qinghai–Tibet Plateau reached a certain height and the inland climate of the Asian continent became arid. In highly differentiated genomic regions, the relative divergence (F_ST) and absolute divergence (d_xy) were significantly higher than the genomic background, θπ and shared polymorphisms decreased whereas fixed differences increased, which indicated that natural selection played a key role in the reproductive isolation of the two species. In addition, we found several genes that were related to reproduction that may be involved in explaining the reproductive isolation. Using phylogenetic trees resolved from haplotype data of Populus tomentosa and P. adenopoda, the maternal origin of P. tomentosa from P. adenopoda was likely to be located in Hubei and Chongqing Provinces.